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Abstract. We employed recently computed evolutionary 
white-dwarf models with helium cores, supplemented by 
heavier models with carbon-oxygen cores, in order to in- 
vestigate the ages of millisecond pulsar systems based on 
the cooling properties of the compact companions. Con- 
trary to the behaviour of more massive white dwarfs, 
the evolutionary speed of low-mass white-dwarf models 
is substantially slowed down by ongoing hydrogen burn- 
ing. By comparing the cooling ages of these models with 
the spin-down ages of the pulsars for those systems for 
which reasonable information about the compact compan- 
ions is available, we found good correspondence between 
both ages. Based on these models any revisions concern- 
ing the temporal evolution of millisecond pulsars do not 
appear to be necessary. 
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1. Introduction 

Millisecond pulsars are thought to be components of low- 
mass binary systems in their final stage of evolution: the 
neutron star which has been spun up by accretion of 
matter from a low-mass evolved companion is now be- 
ing slowed down by emission of magnetic dipole radia- 
tion (recycled radio pulsar). The companion, after having 
transferred most of its envelope mass towards the neutron 
star, remains as a white dwarf of rather a low mass whose 
core consists, in the majority of the known cases, of he- 
lium. The characteristic age, or so-called spin-down age, of 
the (recycled) pulsar depends on the physics how the neu- 
tron star's rotational energy is converted into non-thermal 
emission of electromagnetic energy. On the other hand, 
the white-dwarf age is ruled by the white dwarf's thermo- 
mechanical structure and the transformation of gravother- 
mal energy content into thermal emission of photons from 
the surface. Any age determinations of the pulsar and the 



dwarf component should give the same answer, provided 
our physical understanding of the pulsar's slow-down pro- 
cesses and the white dwarf's cooling properties is correct. 

So far, no general consensus on this matter has been 
achieved. Under the assumption that the cooling proper- 
ties of low-mass white dwarfs are ruled by rather simple 
laws as is known from evolutionary calculations of more 
massive white dwarfs with carbon-oxygen cores (cf. Iben 

1986j Blocker 



& Tutukov 1984; Koester & Schonberner 



1995), large age differences between the pulsars and their 
dwarf companions have been found. In general, the white 
dwarfs appear to be much younger than the pulsars (cf. 
Hansen & Phinney 1998b for a recent, detailed account). 
The best-studied example is the PSR J1012+5307 system, 



for which Lorimer et al. (1995) determined 7 Gyr for the 
spin-down age of the pulsar, but only about 0.3 Gyr for 
the white dwarf's age. Note that the usual spin-down age 
determinations are based on the assumption that the ini- 
tial rotational period after completion of the spin-up by 
accretion is much smaller then the present one, and that 
the pulsar emits magnetic dipole radiation (braking index 
n = 3). A summary of the assumptions inherent in the 
derivation of characteristic or spin- down ages of pulsars is 
given in Hansen & Phinney ( 1998b ). A discrepant result as 
found for PSR J1012+5307, if true, would have important 
consequences for the details of the accretion pro cess a nd 
the following spin-down phase (cf. Burderi et al. 1996). 

A larger sample of millisecond pulsar systems with 
white-dwarf companions has recently been investigated 
by Hansen & Phinney (1998b), using a grid of low-mass 
white-dwarf sequences especially computed for this pur- 
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purse (Hansen & Phinney 1998a). In most cases spin-down 
and cooling ages appeared to be discrepant to various de- 
grees, and the authors were able to constrain the initial 
spin periods and spin- up histories for individual systems, 
especially also for the PSR J1012+5307 system. How- 
ever, the white-dwarf models which this study is based 
on, are generated from ad-hoc assumed initial configura- 
tions. These configurations appear not to be consistent 
with respect to the thermo-mechanical structures and un- 
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processed, hydrogen-rich envelopes with what would be 
adequate for companions in these pulsar binary systems. 

The early investigations concerning the evolution of 
helium white dwarfs made by Webbink (1975) indicated 
that the final cooling is slowed down considerably by on- 
going hydrogen burning via the pp cycle. Obviously the 
cooling behaviour of low-mass white dwarfs depends on 
the size of the still unprocessed hydrogen-rich envelope, 
i.e. whether this envelope is massive enough as to sus- 
tain burning temper atures at its bottom for a long time 
span. The Webbink ( 1975 ) white-dwarf models are, how- 
ever, just evolved main sequence stars without any con- 
sideration of mass loss. 

Since white-dwarf envelope masses cannot be guessed 
from first principles, they must rather be determined by 
detailed evolutionary calcul ations . A step in this dir ection 
was made by Alberts et al. ( 1996 ) and Sarna et al. ( 1998 ) 
who modelled the PSR J1012+5307 system and in partic- 
ular the evolution of the mass giving companion. It turned 
out that the donor shrinks below its Roche lobe while still 
having a rather massive hydrogen-rich envelope which is 
able to keep hydrogen burning dominant even through the 
white-dwarf cooling phase. The evolution was slowed down 
to such an extent that the discrepancy with the spin-down 
age of the pulsar vanished completely. 

Strictly speaking the strength of hydrogen burning, 
and hence the cooling age of an observed white dwarf, 
depends on the size of the envelope before entering the 
cooling path. This envelope mass can be reduced because 
of thermal instabilities of the burning shell when the CNO 
rate dies out, namely by 



— enhanced hydrogen consumption during the instability 
(flash) itself, and by 

— a possible Roche-lobe overflow driven by the rapid en- 
velope expansion. 

The latter case was dominant for the evolution of the 
Iben & Tutukov ( p.986[ ) 0.3 Mq helium white-dwarf model: 
Roche-lobe overflow due to the flash-driven envelope 
expansions reduced the envelope mass below the crit- 
ical value necessary for hydrogen burning. The white- 



— the anti-correlation between core mass and size of en- 



dwarf models of Webbink ( |1975| ) and Sarna et al. ( |1998[ ) 
experienced phases of unstable hydro gen burning for 
M < 0.2 M Q (but see Driebe et al. |l999j for a discussion). 

Recently Driebe et al. ( |1998[ ) published a grid of evolu- 
tionary tracks for helium white-dwarf models which were 
generated by enhanced mass loss applied at different posi- 
tions along the red-gi ant br anch of a 1 Mq sequence (see 



also Iben & Tutukov |l986j , Castellani et al. 1994). This 



method mimicks to some extent the mass transfer in bi- 
nary systems and allows to get reliable post-red-giant con- 
figurations which are very useful for the interpretation of 
observations. Driebe et al. ( 1998 ) covered the whole mass 
range of interest, and they demonstrated that 



velope (cf. Blocker et al. 1997) determines later the 
nuclear activity along the cooling branch, and that 
— thermal instabilities of the hydrogen-burning shell ap- 
pear to be restricted to the mass range of approxi- 
mately 0.2 to 0.3 Mq. 

The absence of thermal flashes below M = 0.2 Mq agrees 
well with the results of Alberts et al. (1996) but disagrees 
with those of Sarna et al. (1998). Nevertheless, the cooling 
times of our models are in excellent agreement with both 
studies. From the given parameters of the white-dwarf 
component in the PSR J1012+5307 system, Driebe et al. 



(1998) determined then its age to be of 6± 1 Gyr, in good 
agreement with the pulsar's spin-down age of 7.0±1.4 Gyr 



(Lorimer et al. 1995). 

The latest effort in a better understanding of the com- 
bined pulsar-white dwarf systems is that of Burderi et 
al. (199S). They took the pulsar spin-down ages at their 
face value and concluded that the standard assumption 
for the white-dwarf cooling (i.e. without nuclear burn- 
ing) complies with the observations, except for masses 
below approx. 0.2 Mq. There are, however, some facts 
that we would like to point out: Burderi et al. (1998) used 



data 'renormalized' to a standard luminosity of 10~ 2 Lq, 
whereby it remains unclear how ages can be renormalized 
if the temporal evolution of the systems is not known a pri- 
ori. Furthermore, they extrapolated existing white-dwarf 
cooling models into mass regimes where they are not valid 
anymore. 

Because of its importance we felt the necessity to re- 
consider the whole issue by utilizing more realistic evolu- 
tionary models for low-mass white dwarfs. We will show 
in the next section that with such models a consistent de- 
scription of those millisecond pulsar binary systems can be 
achieved for which sufficiently accurate data is available. 

2. Pulsar characteristic times and cooling ages of 
their white-dwarf companions 

We started with the sam ple o f millisecond pulsar systems 
used_by Burderi et al. ( |l998 , 



see their Table 1 and our 
Fig. [l]), but made a few changes: the companion mass for 
PSR J 1012+5307 was updated according to Driebe et al. 
( |1998| ), and the systems PSR J1640+2224 and PSR J0437- 
4715 were omitted because of too uncertain pulsar ages. 
For convenience, the relevant data are collected in Table 
0, and all the listed systems are shown in Fig. [l] where the 
characteristic ages of the pulsars are plotted against the 
possible mass ranges of their (white-dwarf) companions. 

The last column in Table |] gives the white-dwarf 
masses according to the binary evolution calculations of 
Tauris & Savonije ( |l999|) . Within these binary calculations 
a relation between the system's orbital period P or b an d the 
white-dwarf mass can be derived (see e. g. Savonije 1987 
and Rappaport et al. 1995). For the systems discussed 
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Table 1. Mass estimates for the white dwarfs and char- 
acteristic ages r c of the pulsars. References are listed in 
the fourth column. The fifth column gives the white-dwarf 
mass from the P or b — -Mwd relation of Tauris & Savonije 



PSR 


Mass limits 
[Mq] 


T c 

[Gyr] 


Ref. 


M W S D " 
M Q ] 


J1012+5307 


0.15-0.19 


7.0 ± 1.4 


1 


0.19 


B0820+02 


0.20-0.64 


0.11 ±0.02 


1 


0.50 


B1855+09 


0.24-0.29 


5.0 ± 1.0 


1,2 


0.26 


J0034-0534 


0.15-0.54 


6.8 ±2.4 


3,4 


0.21 


J1713+0747 


0.27-0.40 


9.2 ±0.4 


3,5 


0.33 



(1998 



Burderi et al. 
van Kerkwijk et al. (|2000|) 
Hansen & Phinney ( L998bj ) 
van Kerkwijk (priv . comm.) 
Camilo et al. ( 199^ ) 



Table 2. Estimates of effective temperature and surface 
gravity for the white-dwarfs in the MSP systems from Ta- 
bic |l|. The fourth column gives the core composition of the 
white dwarf. 



PSR 


Tcff 


[K] 


logs 


core 


Ref. 


J1012+5307 


8 550 


±25 


6.75 ± 0.07 


He 


1 




8 670 


±300 


6.34 ±0.20 


He 


2 


B0820±02 


20 000. 


. 22 000 


6.0. . .7.4 


He 


3 




15 000 . 


.18 000 


7.75 ... 8.0 


C/O 


3 


B1855±09 


7000. 


.9 000 


7.2 ± 0.1 


He 


3 


J0034-0534 


5 000. 


.8 500 


7.2 ±0.5 


He 


3 




^ 4000 


w 7.8 


C/O 


3 


J1713±0747 


4 000. 


.4 500 


7.35 ± 0.05 


He 


3 



1: van Kerkwijk et al. (1996 ;) 
2: Callanan et al. ( |l998 ) 
3: present work 



here the P or b — A/wd relation of Tauris & Savonije (1999) 
predicts masses which are well within the estimated mass 
limits (see Table [l]) . 

It should be emphasized that the characteristic ages 
given in Table [j] and plotted in Fig. [j] are based on certain 
assumptions (see Introduction) which may not be fulfilled 
in all cases. The corresponding systematic errors are dif- 
ficult to assess and cannot be accounted for in this study. 

Also shown in Fig.[l] are (post-red giant) ages of helium 
white-dwarf models taken from Driebe et al. (1998), sup- 



plemented by ages from evolutionar y whit e-dwarf models 
with carbon-oxygen cores (Blocker 1995 ). The ages are 
given for four effective temperatures as to simplify the 
comparisons with the observed sytems: the range between 
4 000 and 20 000 K embraces roughly the estimated effec- 
tive temperatures of the companions of the systems PSR 
J0034-0534, PSR J1713±0747, PSR J1012± 5307 and PSR 
B0820±02 (cf. Hansen & Phinney |l998b| ). Only for the 
PSR J1012±5307 companion exists a rather accurately 
determined effective temperature (8 600 K, van Kerkwijk 

|1998| ). 



et al. 1996, Callanan et al. 



Note that the model ages are counted from the be- 
ginning of the post red-giant phase, i.e. they include also 
the contraction towards the white-dwarf regime. For white 
dwarfs of low mass this contraction from a giant towards a 
white-dwarf configuration makes up for a significant frac- 
tion of their ages and must be accounted for in younger 
systems like PSR B0820±02. 

The temporal behavior of the models, as shown in 
Fig. [I], is determined by the following facts: 

— The compositional differences between the lighter and 
heavier white dwarfs causes the obvious age break 
around Mwd = 0.5 Mq. 

— Hydrogen burning in the helium white dwarfs is, for a 
given temperature, responsible for the strong increase 
in age with decreasing mass. 



— At high temperatures, the lighter models (Ss 0.23 Mq) 
are still in a pre white-dwarf phase very close to the 
turn-around point with rather low (post red-giant) 
ages. 

Fig.|l] clearly demonstrates that for an effective tem- 
perature of 4 000 K our low- mass white-dwarf models ex- 
ceed ages of 10 Gyr, roughly consistent with the pulsar 
characteristic ages. From the positions of individual pul- 
sars (with error bars) we can estimate effective tempera- 
ture and gravity ranges, and also the internal composition, 
to be expected for the white-dwarf companions. The re- 
sults are collected in Table |[ Since for some of the white 
dwarfs temperature estimates based on photometry are 
available (Hansen & Phinney |l998b| ), consistency checks 
are possible. 

2.1. PSR J1012+5307 

The white dwarf in PSR J1012±5307 is so far the best 
studied companion of all known systems (van Kerkwijk et 
al. 1996; Callanan et al. 199S). With its surface parameters 
known, and together with our evolutionary white-dwarf 
models, a consistent description of the whole system is 
found (Fig.g, upper two panels). Since the mass ratio of 
both components is known, the system's inclination can 
be determined as well. The course of the inclination angle 
with white-dwarf companion mass is shown in the lower 
panel of Fig.|| for two limiting pulsar masses. The incli- 
nation can be expected to lie between 40 and 50 degrees. 

2.2. PSR B 0820+ 02 

The mass of the companion is rather ill-defined, and 
it could have a helium or carbon/oxygen core (cf. Ta- 
ble ||) . However, the temperature estimate, 15 250 ± 250 K 
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Fig. 1. Time-scales vs. white-dwarf masses. Boxes indicate millisecond pulsar systems with reasonably known spin- 
down ages and companion masses. The sizes of the boxes indicate the errors listed in Table 1. The lines with symbols 
(+: helium white dwarfs, x: carbon/oxygen white dwarfs) mark the post red-giant ages of white-dwarf models at four 
different temperatures along their cooling tracks taken from Driebe et al. (1998) and Blocker (1995). For more details 
see text. 



(Hansen & Phinney 1998b), allows only a C/O white 
dwarf with at least ~ 0.5 M Q and a rather high grav- 
ity (log g ~ 8), in agreement with the results of Tauris & 
Savonije (]199S|) (cf. Table |). 



2.3. PSR B 1855+09 

The temperature of the companion has recently been pho- 
tometrically determined by van Kerkwijk et al. (2000) to 
be T e s = 4 800+800 K. Given its accurately known mass of 

+0.028 M 
-0.016 M G 

sar timing (Kaspi et al. 



0.258to.oi6 due to t he me asured Shapiro delay of pul 



199J), this low temperature cor- 



responds to a cooling age of 10 Gyr using our model s and 
3 Gyr using the models of Hansen & Phinney ( 1998a ) with 
their smaller envelope masses (< 3 x 10 -4 M Q ). The char- 
acteristic age of the pulsar is, however, 5 Gyr (cf. Table 
|l|) , a fact which on one hand might be a hint for a smaller 
braking index of the pulsar as discussed by van Kerkwijk 
et al. (2000). On the other hand this result strengthens 



the dependence of age on the thickness of the hydrogen 
envelope. 

According to our 0.259 M Q model, hydrogen burning 
ceases at about 10 Gyr, leaving a final unprocessed enve- 
lope of 5 • 10~ 4 M Q . The envelope mass right after the 
last shell flash is ss 2 • 10~ 3 M Q , a value which is, however, 
subject to uncertainties like flash strength and metallic- 
ity. In additional calculations we artificially reduced this 
envelope mass by invoking mass loss on the upper cool- 
ing branch and followed the evolution of the models in 
the usual manner. It turned out that an earlier reduction 
of the envelope mass to ~ 5 ■ 10 _4 M Q after m 0.5 Gyr 
(at T c ff w 10 4 K) would be sufficient to give a cooling age 
of 5 Gyr at the desired effective temperature of 4800 K. 
We note that with this reduced envelope mass, hydrogen 
burning becomes insignificant below T c g ph 10 4 K. 



The mass of the chemically homogeneous envelope is de- 
fined as the total mass of the layers above the hydrogen ex- 
hausted core. 
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J1 01 2+5307: 7.0 +/- 1.4 Gyr 




0.55 



0.3 0.35 0.4 

M WD/ M sun 



0.55 




0.55 



M W n/M =l 



Fig. 2. Upper panel: Isochrones for three white-dwarf 
ages, 7. Oil. 4 Gyr, in an effective temperature vs. mass di- 
agram. The effective temperature of the white-dwarf com- 
panion of PSR J1012+5307 is indicated by the horizon- 
tal dashed lines (long dashed: van Kerkwijk et al. 1996; 



short dashed: Callanan et al. 1998). Middle panel: The 

same isochrones but in a surface gravity vs. mass diagram. 
Again the observed values for PSR J1012+5307 are also 
given (dashed lines). Lower panel: The PSR J1012+5307 
system's inclination for two pulsar masses vs. the compan- 
ion mass. 



2.4. PSR J0034-0534 



Hansen & Phinney (1998b) estimated a very low temper- 
ature limit of T e g < 3 500 K for the white dwarf, yielding 
ages of more than 10 Gyr if helium models are used, which 
is to be compared with the pulsar's characteristic age of 
6.8 ± 2.4 Gyr (Fig. |l|) . Consistency between the pulsar's 
and the white dwarf's age can only be achieved if we as- 
sume the white dwarf to have a carbon/oxygen core and 
a mass of » 0.5 Mq. Such a model cools considerably 
faster and reaches 3 500 K well within about 6 Gyr. We 
note that this mass is noticeably larger than estimates 
of other studies. The P or b — Mwd relation of Tauris & 



Savonije (1999) gives M ta 0.21 M , and the photometric 
measurements of Lundgren et al. ( |l996a| ) M « 0.23 Mq. 
The cooling models of Hansen & Phinney ( |1998b ) pre- 
dict Mwd ~ 0.32 Mq as an upper limit. The discrepancy 
between these results and our mass estimate (C/O- white 
dwarf) might be related to the same problem as in the case 
of PSR B1855+09. If the companion is a helium-white 
dwarf with M 0.2 . . . 0.3 M it is prone to hydrogen 
shell flashes with the corresponding uncertainties. 

2.5. PSR 1713+0747 



Hansen & Phinney ( |l998b[ ) give a white-dwarf effective 
temperature of 3 400 ± 300 K which is just slightly lower 
than the one predicted by our helium models (see Ta- 
ble ||) . The white-dwarf mass is not very sensitive to the 
temperature value and consistent with the es timate of 
Mwd ~ 0.33 M Q from Tauris & Savonije fll999|). 



2.6. Other millisecond pulsar systems 

In addition to the sample discussed above we investigated 
some other MSP systems with respect to the possible 
(g, T c ff ) combinations for the white-dwarf components (see 
Table ||) . The results are illustrated in Fig. || where ef- 
fective temperature and surface gravity are plotted as a 
function of the white-dwarf mass for the systems listed in 
Table ||. The white-dwarf masses can be estimated from 
the P rb — Mwd relation of Tauris & Savonije (199E). Tak- 
ing these masses, it is straightforward to determine effec- 
tive temperatures and surface gravities of the white-dwarf 
companions (see Table ||). 

We find the systems in the first row of Fig. |^ 
(J1640+2224, J0751+1807, J1045-4509) to be consistent 
with He- WD companions if T off > 4 000 K. At this tem- 
perature our helium models reach cooling ages comparable 
with the age of the galactic disk (w 10 Gyr). For the white 
dwarf in J1640+2224 a temperature estimate is available: 
4 500 ± 1100 K (Hansen & Phinney |1998b| ). This value im- 
plies, for 10 Gyr, a white-dwarf mass of ps 0.35 ± 0.05 Mq 
with a surface gravity of log g « 7.6 (cf. Fig.||). This mass 
i s in g ood agreement with the result of Tauris & Savonije 
(fHiJ (MTS99 «0.37M o ). 

For the other systems the temperature ambiguity 
(more than one white dwarf mass for a given tempera- 
ture, see. Fig. |j) does not allow to exclude the compan- 
ions to be C /O-wh ite dwarfs, but the mass values of Tauris 
& Savonije ( 1999] ) indicate that all these systems should 
contain helium white dwarfs. 



3. Conclusions 

From the present paper, together with previous efforts 
which concentrated solel y on the PS R J10 12+5307 sys- 
tem (Alberts et al. 1996 ; Sarna et al. 1998 ; Driebe et al. 
1998 ), it becomes obvious that a consistent description 
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of the millisecond binary systems with compact compan- 
ions can only be achieved by using evolutionary model 
calculations of white dwarfs which include their complete 
pre- white-dwarf history. The key to the solution of the ap- 
parent age paradoxon between the pulsar and its white- 
dwarf companion is the fact that low-mass white dwarfs 
have massive, still unburnt envelopes that sustain hydro- 
gen burning at their bases for a long time. Hydrogen burn- 
ing slows down the cooling of a low-mass white dwarf to 
such an extent that cooling ages become comparable to, 
or may even exceed, observed pulsar spin-down ages. 

Employing our evolutionary helium white-dwarf mod- 
els, supplemented by those with carbon-oxygen cores, 
we demonstrated that, next to the already studied PSR 
J1012+5307 system, also in other millisecond pulsar bi- 
nary systems with reliable information on pulsar age and 
companion properties, as mass and temperature, reason- 
able agreement between the components' ages is achieved. 
The use of white-dwarf models with ad hoc assumed en- 
velope masses may lead to erroneous interpretations since 
these envelope masses are usually much smaller than those 
which follow from complete evolutionary calculations. 

It appears to us that upon using realistic white-dwarf 
models in interpreting millisecond binary systems there is 
no need to modify existing ideas of the spin-down process 
of pulsars. Clearly a larger sample of well-studied systems 
like PSR J1012+5307 would be very important in investi- 
gating more precisely the cooling theory of white dwarfs 
and the braking of radio pulsars. 

Acknowledgements. We would like to thank Marten van Kerk- 
wijk and Gerrit Savonije for helpful discussions and comments. 
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Table 3. Data for selected millisecond pulsars. P S pin: pulsar's spin period; P: change of Pspin! t c : characteristic age; 
f(M): mass function; M^j^ 9 : white-dwarf mass according to the P or b — Mwd of Tauris & Savonije ( 1999| ); T Q g,g: 
predicted effective temperature and surface gravity for M\vd = Mwd 9 (compare Fig. ||) . 
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1: Taylor et al. (1993 , 199 5 (unpubl. work, available via anonymous ftp from pulsar.princeton.edu (128.112.84.73))) 
2: Lundgren et al (ll996b| ) 
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4: Sandhu et al. (jpL997|) 
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Fig. 3. Possible effective temperatures (left y-axis, solid lines) and surface gravities (right y-axis, dashed lines) as a 
function of white-dwarf mass for the millisecond pulsar systems of Table || for three different cooling ages. The middle 
lines belong to the given characteristic age (see Table ||), except for J1640+2224 where the white-dwarf properties are 
shown for r = 10 Gyr. The upper and lower lines refer to ages r c ± 2 Gyr. For J0437-4715, where only an upper limit 
for t c is known, lines are given for r = 1, 3 and 6 Gyr. Note that the lower age limit is represented by the upper lines 
for T e s and the lower lines for g. 



